Theory of phase fluctuating d-wave superconductors and the spin response in 

underdoped cuprates 
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The minimal theory of spin of gapless quasiparticles coupled to fluctuating vortex defects in 
the phase of the d-wave superconducting order parameter at T = is studied. We find a single 
superconductor-spin density wave phase transition, which may be fluctuation induced first-order, at 
which vortices condense and the chiral symmetry for fermions dynamically breaks. We compute the 
spin-spin correlation function in the fluctuating superconducting state, and discuss some prominent 
trends in the neutron scattering data on underdoped cuprates in light of our results. 



The underdoped high temperature superconductors 
are highly anisotropic, quasi two-dimensional materials, 
in which thermal and quantum fluctuations of the phase 
of the superconducting order parameter (OP) should be 
important well below the pseudogap temperature 0. It 
has recently been shown that a phase fluctuating d-wave 
superconductor (dSC) at zero temperature (T = 0) is in- 
herently unstable towards the formation of the insulating 
spin density wave (SDW) state with the loss of phase co- 
herence Q, Q. This result follows from the realization 
that the effective low-energy theory for the gapless d- 
wave quasiparticles, besides the usual spatial symmetries, 
also possesses an additional internal ("chiral") symme- 
try. This symmetry becomes dynamically broken when 
the superconducting phase coherence is lost via prolifer- 
ation of vortex defects. The effective theory of the phase 
fluctuating dSC in this formulations is closely related to 
the three dimensional quantum electrodynamics (QED3) 
0, in which the "charge" is proportional to the vortex 
condensate 0, H- In first approximation the fluctua- 
tions of the gauge-field may be completely neglected in 
the superconducting state, which yields quasiparticles as 
well defined excitations. When the vortices condense and 
the charge in the effective QED3 becomes finite, the chi- 
ral symmetry of the dSC becomes dynamically broken 
|{|. In the present context this translates into the SDW 
ordering, with confined spin-1/2 (spinon) excitations Q. 

Several issues of direct relevance for underdoped 
cuprates naturally arise in such a theory of the phase fluc- 
tuating dSC. Can the SC and the SDW long-range orders 
coexist? Several intriguing recent experiments [g , , [3 , 
find possible signs of such coexistence of the two inimical 
orderings. If there is a single dSC-SDW quantum phase 
transition, on the other hand, what should be its univer- 
sality class? What are the effects of the vortex fluctua- 
tions inside the superconducting state? These questions 
all call for a better understanding of the quasiparticle- 
vortex interaction, particularly in the superconducting 
state. In this Letter we investigate probably the simplest 
theory of the low-energy d-wave quasiparticles coupled 
to the complex scalar vortex field (dual to the supercon- 
ducting OP). The theory describes the decoupled spin 
sector of the phase coherent dSC at T = and sub- 
ject to strong quantum vortex fluctuations. In under- 
doped high temperature superconductors such fluctua- 



tions should presumably arise from the Coulomb interac- 
tion that becomes effectively stronger towards half filling, 
and tends to disorder the phase of the superconducting 
OP. If extended to the non-superconducting phase, our 
theory reduces to the previously studied QED3. First, 
we argue that the condensation of vortices (i. e. the loss 
of superconductivity) and the breaking of the chiral sym- 
metry for fermions (the SDW instability) in our theory 
coincide. Second, we show that the quantum dSC-SDW 
phase transition can be either in the modified XY uni- 
versality class, or it could be fluctuation-induced first- 
order, depending on the values of the couplings in the 
theory. Finally and most importantly, we calculate the 
spin dynamics induced by vortex fluctuations deep inside 
the dSC in the leading approximation. The intricate evo- 
lution of the spin response with frequency observed in in- 
elastic neutron scattering experiments on YBCO |fj|.[Tfj|. 
|llj follows quite naturally from our results. In particu- 
lar, we predict the appearance of four weak and narrow 
"diagonally" incommensurate peaks at low energies, the 
energy of which vanishes with the superconducting T c . 

Consider the quantum mechanical (T = 0) action 
for the low energy quasiparticles in the two-dimensional 
phase fluctuating dSC, S = J d 3 x£, x = (t, r), and 

N 

£ = V%(3 )1 -» (1 )* i + -»-(Vxi)+ (1) 

*■ — ' IT 

i=l 

|(V-iI)$| 2 + M 2 |$| 2 + ^| 4 . 

Here the fluctuating complex field $ describes the vor- 
tex loops, and (3>) = implies the superconducting 
phase coherence. Two (N — 2) four-component Dirac 
fermions describe the gapless, neutral, spin-1/2 (spinon) 
excitations near the four nodes of the superconducting 
order parameter, one Dirac field for each pair of diago- 
nally opposed nodes |2| . a is the gauge-field that results 
from properly absorbing the singular part of the super- 
conducting phase due to vortices into the spinon fields 
0, 0- 7ii, H — 0,1,2 are the Dirac gamma matrices 
0- Finally, A is an additional auxiliary (Chern-Simons) 
field which facilitates the statistical spinon-vortex cou- 
pling, as will be explained shortly. The single tuning 
parameter /1 2 may be understood as related to doping, 
/1 2 oc (x — x c ), with x c being the critical doping in the 
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underdoped regime. The coupling b describes the short- 
range repulsion between the vortex loops. In our units 
h = c = vp = va = &b = e = 1- 

We may motivate the field theory (1) as follows _12j: 
the integration over the gauge-field a leaves the spinon- 
vortex coupling iA ■ J$, where V x A = irJ^, and 
and J$ are the spinon and vortex current densities, re- 
spectively. Vortices and spinons therefore see each other 
as sources of magnetic flux, and circling around a vor- 
tex with a spinon (or vice versa) leads to a phase change 
of 7T. This is precisely the statistical (Aharonov-Bohm) 
part of the interaction between vortices and spinons , 
13 1 E2 ■ Without fermions (N — 0) , the integration over 
the gauge fields leads to the standard |$| 4 theory, in the 
universality class of the XY model. Strictly speaking, for 
N = the theory should reduce to the Higgs scalar elec- 
trodynamics, dual to the |$| 4 theory. This can be accom- 
plished by introducing an additional fluctuating gauge 
field into the theory that would mediate a long-range 
interaction between vortices and represent the electrical 
charge. Following the previous work 0, H, 0] we will 
assume spin and charge to be separated at low energies, 
so we may neglect the charge degrees of freedom. This 
assumption is surely justified inside the superconducting 
state which is the subject of the present work. For subtle 
issues that arise in the normal state the reader is referred 
to (l4|. Finally, the velocity anisotropy (vf ^ va) and 
all the local interactions between quasiparticles |2| have 
been neglected, as irrelevant at low energies |l5j. |16||. 

Let us discuss first how the theory (1) is related to the 
QED 3 @, 3, in the simplest mean-field approximation 
for the vortex field. Neglecting the fluctuations in $, 
in the superconducting phase ({$) = 0), the integration 
over A constrains V x a = 0, and spinons are essen- 
tially free. When vortices condense (($) 0), on the 
other hand, the gauge-field A acquires a mass via Higgs 
mechanism, |(<&)| 2 ^4 2 . The Gaussian integration over A 
produces then the Maxwell term for a,~(Vx a) 2 /|($}| 2 . 
Together with the Dirac Lagrangian this constitutes the 
QED3 for spinons, in which the fermions' chiral sym- 
metry generated by 73 and 75 is dynamically broken by 
the generation of the mass term ~ M^^i |f|. In the 
present context such a mass M ~ (^^i) ~ K^)! 2 is pro- 
portional to the SDW OP, with the (incommensurate) 
ordering wave vectors that connect the diagonally op- 
posed nodes of the superconducting OP 2j. Carefully 
including the charge degrees of freedom in (1) one also 
finds the SDW to be an electrical insulator [l4| . 

In the mean-field approximation, the vortex condensa- 
tion and the SDW transition in the theory (1) coincide. 
It is not obvious that this feature survives the inclusion 
of the fluctuations, and it is conceivable that the SDW 
transition may occur within the superconducting phase 
0> E3- To examine this issue we first ask what would 
be the exact propagator for the gauge field a if there were 
no fermions, i. e. when TV = in (1). We then approxi- 
mate the interacting action for a that would result from 



the integrations over $ and A with the effective Gaus- 
sian term that reproduces that exact propagator. In the 
Landau gauge such a propagator is 



Gla.uviPi m ) 



7T 2 U AA {p,m) 



(2) 



where Haa(p, to)(<5 m „ — p^p v ) is the transverse current- 
current correlation function in the l^ 4 theory. Here 
to 2 = /i 2 + 0(b) is the fully renormalizcd "mass" of the 
vortex field. In general, Haa{p, to) = npF + (m/p), for 
to 2 > (superconductor). To the lowest order in b 

F + (z) = — [(4z 2 + l)arctan( — )-2z] + 0(b). (3) 
8n 2z 

Here we generalized our model to the one with n complex 
vortex fields, n = 1 being the case of physical interest. 

We are now in position to study the chiral symmetry 
breaking for fermions with (2) serving as the bare (with- 
out fermion polarization) gauge-field propagator in the 
QED3. Consider the standard large-N Dyson equation 
[3 for the fermion self-energy 



d 3 p Ggg^ip - q,m)Y,{p) 

(27r) 37At V " 1 ' 



P 2 + S(p) 2 



Here [G aa ^{p,m)} 1 = Hf 0)|U/ (p) + [G° aa ^ v {p, m)] \ 
with the one-loop fermion polarization n^ a „„ (p) = 
(Np/8)(S IJ , U — pupv). It is useful to consider first the 
point of the superconducting phase transition m = 0, 
where Eq. (3) implies n^^O) = up(l/16 + 0{b c )), 
where b c is the fixed point value of b in d = 3 [Tsj . Linear 
dependence of n^^p, 0) on momentum is an exact result 
[l9l | . Inserting this into the Dyson equation we find that 
at to = the effect of vortices is only to increase the co- 
efficient N — > N e f f = N + 128/(7r 2 n) + 0(b c ). Recalling 
that the non-trivial solution of the Dyson equation (4) 
exists only for N eff < N c = 32/vr 2 5], for TV = 2 we find 
that chiral symmetry at m = would be already broken 
only for n > n c — 10.44(1 + 32/(97r 2 n)), where we have 
also included the known Q(b r ) correction to Haa(p, 0) in 
the large-n approximation |l8j . Since n = 1 in the phys- 
ical case we conclude that right at the superconducting 
critical point SDW order is most likely absent. 

Little further thought shows that the above result im- 
plies that the SDW order is absent for all to 2 > 0. Indeed, 
this is to be expected, since G° Q (p, to) < G„ n (p, 0), and 
a is only stiffcr in the superconducting phase. To prove 
this, assume that for m 2 > a non-trivial solution of the 
Eq. (4), does exist. Such a solution would then 

satisfy 



t(q) < -Tr 



d 3 p Gaa^uip ~ £?,0)E(p) 



(27T) 



p 2 + S(p) 2 



(5) 



On the other hand, we already established that there is 
only a trivial solution of the Dyson equation at to = 0. 
This means that assuming a candidate function X(p) and 
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inserting it under the integral in the Dyson equation (4) 
for m = will produce only a smaller new £(<?), since un- 
der iterations the physically acceptable self-energy must 
approach the trivial solution. This further implies that 
any E(q) from the domain of attraction of the trivial so- 
lution at m = has to satisfy the inequality opposite to 
(5). A non-trivial £(<z) at m 2 > therefore can not exist 
if it did not exist at to 2 = 0. 

For m 2 < superconductivity is lost, and Haa(p, to) = 
npF-(\(<f>)\ 2 /p), with F-(z) -> 1/16 for z < 1, and 
F-(z) = 2z, for z ^> 1 (Higgs mechanism), to the lead- 
ing order. Since G° aa (p,m) > G® a (p, 0) for to 2 < 0, the 
above argument no longer applies. In fact, it is easy to 
see that there is immediately a non-trivial solution of the 
Eq. (4) when to 2 < 0. Since E(g) > only for q ~ |(4>)| 2 
and smaller, |($)| 2 serves as the effective ultraviolet cut- 
off in the Eq. (4). The Dyson equation then reduces to 
the standard one in the QED 3 [5j. For N = 2 < N c , 
S(g = 0) ~ K^)! 2 , and the chiral symmetry is dynami- 
cally broken. There is no intermediate (quantum disor- 
dered) phase in between the SDW and the dSC in the 
theory (1), unless the exact value of N c in the QED3 is 
actually less than two po|. 

The above argument strongly suggests that there is a 
single dSC-SDW transition in the theory (1), but does 
not say anything about its nature. To study this issue it 
is better to proceed in the opposite direction, and inte- 
grate over the fermions first. This gives dynamics to a via 
the fermion polarization bubble. If we neglect the quar- 
tic and the higher order terms and perform the Gaussian 
integral over a, the result is the Maxwell-like term for 
A; {A ll (p)A u (-p)) = (NTr 2 /8\p\)(5 fJ , v -p ll p u ). Note that 
since the inverse of the above av erag e is non-analytic at 
p = 0, it can not renormalize |2l|. and therefore the 
number of fermions N represents an exactly marginal 
coupling. Assuming a constant <E> in (1) we may further 
integrate A to find the energy per unit volume to be 



07T 



ir 2 N 



(6) 



; 384 v 



2 A$ 2 



), 



where we have rescaled $ 2 /A — > $ 2 , with A being the 
ultraviolet cuttof. Standard analysis of the Eq. (6) 
shows that there is a discontinuous transition for N > 
(4/7r)y26/A). Using the lowest order fixed point value 
b/A = (27r 2 /5) (4 - d) + 0((4 - a 1 ) 2 ) in d = 3 as a crude 
estimate of this bound, we find the first-order transition 
for N > 3.58, in rough agreement with |12j . Of course, 
this conclusion is to be trusted only for W > 1, when 
the first-order transition occurs at large fi 2 , at which our 
neglect of fluctuations in <f> in arriving at the Eq. (6) 
becomes justified. For N <C 1, on the other hand, the 
phase transition in the theory (1) should be continuous, 
with weakly modified XY exponents: v = v xy + O(N), 
V = Vxy + O(N) [lj. The situation is reminiscent of the 
Ginzburg-Landau superconductor [19j , with N 1 anal- 



ogous to the strongly type-I, and N <C 1 to the extreme 
type-II case. 

Finally, we turn to spin dynamics deep inside the dSC, 
when to/6 3> 1. For small momenta, p -C m, we find 
then that G° aa ^(p,m) = ((tt /2±m) + 0(p 2 )){5 tiv -p tl p v ). 
Integrating out such a massive a leads to the effective 2+1 
dimensional Thirring model for fermions in the dSC 



72to 



(7) 



with the summation over repeated indices assumed. Us- 
ing the Hubbard-Stratonovich transformation we can 
rewrite this as 



18m 



Tr[M^M^\ + Tr[M^^M (8) 



Making an ansatz M % J (x) = (M(x)/3)Sijj fi at the saddle 
point, and then integrating out the fermions yields 

/94m 
d 3 x[—^M 2 (x) - Tr[ln(j^ - M(x))}}. (9) 

Expanding further in powers of M(x) we finally write 



S 
N 



d 3 q ,24irm — A loL n, . „, . % 

( 5 + m)M 2 (q) + 0(M i ), (10) 



(2tt 



where M 2 (q) = M(q)M(—q). On the other hand, using 
the definition of the Dirac fields in terms of the electron 
creation and annihilation operators 

2 

(M(x)) = _^-(4(f,T))^cos(2^-r), (11) 



where the vectors ±ifj, i = 1,2 denote the positions of 
the four nodes of the d-wave order parameter. (M(q = 
0)) is therefore the static SDW OP @. Similarly, 



<M 2 (g») 



48Ato 



(12) 



£1 

(^ro") 2 E«^(9 + ™ h u)S z {-q- 2K h - w )> 



(S z (q-2K i) u;)S z (-q + 2K i ,-u))). 



Since the dSC is rotationally invariant, the spin-spin cor- 
relation function is diagonal, {S a (k 1 Lu)Sfj(—k, —u>)) = 
x(k,oj)8 a f3- We may therefore finally deduce that the 
imaginary part of the spin response function x(k,to) = 
x'(k,u) + ix"{k,uj) in the phase fluctuating dSC is 

X "(2K i ± q,uj) = (^H) 2 Im(M 2 (q,to)). (13) 

Analytically continuing to real frequencies iui — > u>, in the 
Gaussian approximation to Eq. (10) we finally obtain 



X {2K i ±q,u) = — — — -j-^- — r , (14 

7r (192to/7t)^ + ur — q l 
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where we have taken m A deep inside the dSC. With 
?7i — > oo Eq. (14) reduces correctly to the non-interacting 
limit, in agreement with p^ . 

First, let us connect the uniform vortex susceptibility 
m to the measurable superconducting T c . Assuming a 
continuous dSC-SDW transition, m 2 oc (/i 2 — /i 2 ) 7 , with 
7 = i>(2 — rj) characterizing the quantum critical point; 
also, T c oc (/i 2 — Hc) zl/ ■ Since the dynamical critical ex- 
ponent in the theory (1) is z = 1, m oc Tt in/2) - Sum- 
ciently away from the critical point we may neglect the 
anomalous dimension 77 and find the result dictated by 
the (engineering) dimensional analysis, 192to/?t = cT c , 
where c is a number. 

Next, we look for the position of the maximum of the 
spin response function x"(<Tj w ) a t a fixed energy to. Eq. 
(14) implies that for u <C cT c , the maximum values are 
located at four "diagonally" incommensurate wave vec- 
tors ±2X1.2 (i. e. at q = 0). As frequency increases the 
peak intensity grows, and recalling that vp 3> and 
taking the lattice periodicity into account one finds the 
four peaks overlapping first at four " parallel" incommen- 
surate positions [22J. With the further increase of fre- 
quency, at some point four initial peaks start overlapping 
at (tt, 7r), and for a while the "commensurate" response 
dominates. At lu > cT c the maximum in Eq. (14) shifts 
to a |<f| 7^ 0, which implies a weak redistribution of the 
commensurate peak to four "parallel" incommensurate 
positions at largest frequencies. The predicted evolution 
of the spin response is in q ualitative agreement with the 
observations in YBCO [ll| . More detailed and quantita- 



tive picture will be presented in a future publication. 

Consider then the commensurate response at k = 
(7r,7r), that in our notation corresponds to some q =/= 0. 
As a function of frequency, at T — x" vanishes below 
some cutoff energy lu c — \q\. Also, the maximum response 
is at the energy luq = \/\q\ 2 + (cT r ) 2 , which decreases 
with decreasing T c , as observed 0, [l(j. As T c — > 0, 
however, the commensurate peak energy u>o — > \q\ 0. 
Such a finite luq (« 20meV) as T c — » is in agreement 
with the data in YBCO (Fig. 29 in the second reference 
0). Remarkably, the cutoff energy that can be estimated 
from a different set of data on YBCO (Fig. 17 in [lj| ) 
is also lu c ps 20meV, consistent with our prediction. At 
large energies, x" should behave as ~ 1/co, which also 
appears to be in general accord with the data. 

Finally, at the incommensurate (q — 0) wavevector, 
x" (2Ki,oj) ~ lu at small energies, with the peak at 
loq = cT c . Whereas the energy of the resonance should 
go to a finite value as T c vanishes with underdoping, we 
predict that the energy of the "diagonally" incommen- 
surate peaks, which should become discernible at lower 
energies, should extrapolate to zero. Detection of such a 
soft mode inside the dSC state would provide the smok- 
ing gun evidence for (1) as the low-energy theory for spin 
of underdoped cuprates. 
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